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Magnetization and susceptibility studies on single crystals of the pyrochlore Ho2Ti2 07 are re- 
ported for the first time. Magnetization isotherms are shown to be qualitatively similar to that 
predicted by the nearest neighbor spin-ice model. Below the lock-in temperature, T* ~ 1.97 K, 
magnetization is consistent with the locking of spins along [111] directions in a specific two-spins-in, 
two-spins-out arrangement. Below T* the magnetization for S|[[lll] displays a two step behavior 
signalling the breaking of the ice rules. 
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Materials containing antiferromagnetically coupled 
magnetic moments which reside on geometrical units, 
such as triangles, that inhibit the formation of a collinear 
ordered state, often display phenomena known broadly 
as geometrical frustration [l]]. A pyrochlore lattice 
with corner-sharing tetrahedra of magnetic ions running 
through it, is an ideal system for studying geometri- 
cal frustration. The discovery by Harris et al. D that 
a magnetic pyrochlore with ferromagnetic interactions 
and strong single ion anisotropy could be just as frus- 
trated as the analogous antiferromagnets is counterintu- 
itive The prototypical example of these ferromag- 
netic compounds is Ho 2 Ti 2 07 (|,||||, however this has 
been questioned recently 0. Such materials are termed 
spin-ices because their physics maps directly onto the 
historic problem of proton disorder in hexagonal ice. 

The magnetic ground state of the ferromagnetic py- 
rochlore with Ising anisotropy along the [111] axis con- 
sists of two spins pointing in and two spins pointing out 
of the individual tetrahedra ||. This is identical to the 
lattice formed by the mid-points of the oxygen-oxygen 
bonds in ice This resulted in the "spin- ice" label, 
and is one of two conditions required for this class of 
materials, the other being a zero-point entropy close to 
(1/2) ln(3/2) which has been measured in Dy2Ti2C>7, and 
is identical to the proton ordering in real ice. The discov- 
ery of spin-ice prompted considerable research by other 
groups j7|]9|]To|] , however until recently neither Dy2Ti2C>7 
nor Ho2Ti207 have been shown to meet these two crite- 
ria. Due to the neutron absorbing nature of dysprosium, 
neutron scattering measurements on Dy2Ti2C>7 to inves- 
tigate the nature of the magnetic ordering have yet to be 
performed. Recently, Bramwell et al. reported the first 
low temperature C m results on Ho2Ti 2 07 Q. In these 
data a maximum is seen in the magnetic contribution, 
C m , to the specific heat at 1.97 K which we will refer to 
as the blocking temperature. Their data, heat capacity 
and neutron diffraction, clearly shows that Ho2Ti207 re- 
quires a spin-ice model that includes longer-range, dipo- 
lar interactions when discussing its properties down to 
300 mK. 

To investigate the bulk magnetic properties of 
Ho2Ti 2 07, we have performed both a.c. and d.c. mag- 
netization measurements on single crystals. These are 
the first single crystal measurements, with fields as large 
as 9 T applied along the three principle cubic crystallo- 
graphic axes [100], [110] and [111], of this geometrically 
frustrated materials. These data are in excellent qual- 
itative agreement with calculations J5|; in particular a 
step is seen below the lock-in temperature, T* when the 
magnetic field is applied along the [111] axis. Unlike the 
predictions, there is little difference in the low temper- 
ature saturation magnetization for the different crystal 
orientations. 



The single crystals were synthesized using the floating- 
zone method described elsewhere |ll]]. Polycrystalline 
samples of Ho2Ti207 were prepared by firing, in air, sto- 
ichiometric amounts of high purity (> 99.99%) H02O3 
and TiC>2 for several days with intermittent grindings to 
ensure a complete reaction. The polycrystalline material 
is then used to grow large (> 1 cm 3 ) single crystals in 
an infrared image furnace with the aid of a small seed 
crystal. A small piece of the resulting large single crys- 
tal was ground up for characterization by X-ray powder 
diffraction (XRD). The impurity content is less than 3% 
(limit of our XRD) and the room temperature cubic lat- 
tice parameter is 10.102(2) A. 

Magnetization measurements between 1.8 and 300 K 
were performed in a fully calibrated Quantum Design 
PPMS system inside a 9 T superconducting magnet. The 
zero field susceptibility measurements used a 10 Oe AC 
field at a frequency of 1000 Hz. The high field magne- 
tization measurements were made with a DC extraction 
technique. 

Previous magnetic measurements on Ho2Ti207 have 
been interpreted in different ways |7],[l2],[l3]] . We will dis- 
cuss our results within the framework of a ground-state 
spin one- half doublet S = 1/2 with an effective g value 
of g 1 = g and g± = where parallel is the [111] direction 
corresponding to the direction of the applied field. We 
will use the standard form |14| 



1+ T 



S(S+l)J/3, (1) 



where J is the effective interaction between the Ho spins. 
As noted by others Jll| , 9 is going to be the sum of the 
exchange, dipole, and crystal field contributions. We do 
not attempt to extract these individual contributions to 
9. The results of the AC zero field susceptibility mea- 
surements are shown in Fig. |l|. For T > 10 K, the 
data are consistent with that expected for g — 20 and 
weak ferromagnetic interactions as previous report [ OJ . 
In the range 3-10 K, the data can be described well using 
9 = —1.8 K which corresponds to an antifcrromagnetic 
interaction, in agreement with the results of Siddharthan 
et al. 0. Jana et al. Jl3| and other have already noted 
that the Curie- Weiss temperature is very dependent on 
the range of temperatures studied, this is undoubtedly 
due to crystal field contributions to the effective inter- 
action. Rosenkranz et al. |l6| ] have shown that the low- 
est crystal field excitation is at ~ 200 K. If we fit our 
high temperature, T > 200 K, DC magnetization data, 
(Fig. |), with B||[lll]=l kOe, we find 9 w +1.6 K, in 
agreement with others From these results, we con- 
clude that the Ho moments can indeed be considered as 
5 = 1/2 spins, with g = 20 and weak, but dominant 
ferromagnetic interactions in the paramagnetic state. It 
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should be noted there that this net ferromagnetic inter- 
action between Ho spins is a combination of antiferro- 
magnetic exchange interactions and ferromagnetic dipo- 
lar coupling. 

Although not shown here, it can be shown that C m is 
also well described by the nearest neighbor spin-ice model 
at temperatures above IK. Using the data from Bramwell 
etal. H we have calculated the magnetic entropy AS m 
by integrating C m /T and the results are shown in the 
inset of Fig. |[ The error bars represent the estimated 
uncertainty in the entropy, mainly due to the fit of hyper- 
fine component. The solid line is the entropy expected 
from Harris et al. Q. The behavior is remarkably simi- 
lar to that observed in Dy2Ti2C>7 || where the applica- 
tion of a 1 T magnetic field nearly recovers the missing 
(l/2)i?ln(3/2) of entropy expected for a spin 1/2 system 
(_Rln2). This is in sharp contrast to a previous report 
M , and is consistent with Ho2Ti207 being a spin-ice com- 
pound. 

One naively expects magnetic ordering for T < \9\. 
However, the geometric frustration leads to the absence 
of magnetic order at least down to 0.05 K ||. Though 
long-range magnetic order does not appear, the magnetic 
interactions can lead to short-range local ordering. Mea- 
surements of the DC magnetization as a function of mag- 
netic field are shown in Fig. ||. Before each of these 
measurements, the temperature was raised to 10 K in 
zero applied field before returning to the desired tem- 
perature and applying a magnetic field. For T > 2 K, 
the curves for £?||[111] and B||[100] are identical, while 
for T < 2 K drastically different behavior is observed. 
The data for the S||[110] orientation is rather interest- 
ing. While the other two orientations show identical be- 
havior for T > T* , suggestive of isotropic magnetization, 
the B||[110] data is rather different. Also, for T < T* 
the J5||[110] data does not differ appreciably from the 4 
K data. In other words, the locking of the spins that 
we believe to occur for T < T* seems to have little ef- 
fect on the behavior when _B||[110]. It should be noted 
that when B||[110], two of the four spins are perpendic- 
ular to the field and cannot order. In the other direc- 
tions all spins have a component of their moment in the 
field direction. In general the data can be explained by 
the spins having strong Ising anisotropy along the [111] 
axis. Below T* ~ 1.9 \6\ the spins lock into a spe- 
cific set of [111] direction so that a tetrahedron has the 
two-spins-in, two-spins-out configuration. From Monte 
Carlo calculations of the nearest neighbor spin-ice lat- 
tice H, it was shown that the degeneracy breaking oc- 
curs via a two step process for B||[lll]. First the spin 
pointing along the [111] direction is aligned with the field 
leading to an average moment in the field direction of 
(l+l/3+l/3-l/3)/4=l/3 followed by a breaking of the 
two-spins-in two-spins-out ice rules with an average mo- 
ment of (l+l/3+l/3+l/3)/4=l/2 at higher fields. This 
predicted behavior yields a plateau at 3.33 /ib/Ho atom 
followed by saturation at 5.00 /j,b/Ho atom for the case 
g = 20. Our data for £?||[111] clearly shows a plateau at 



~ 3.3 /ie/Ho appearing as temperature is lowered, which 
agrees well with the predicted behavior. However, the 
saturation magnetization is approximately 20% higher 
than calculated. For B||[100] and 5||[H0], the data at 
T = 1.8 K qualitatively agree well with the predicted be- 
havior In particular, for T < T* , the field required 
to reach saturation is significantly less for £?||[100] as for 
5||[110]. It is seen that, in stark contrast to the pre- 
dicted behavior of a saturation moment of 5.78 /ib/Ho 
for B||[100] and 4.08 ^ B /Ho for B||[110], the measured 
saturation moment of 5.9 /ib/Ho is close to the value 
expected for £?||[100] and is nearly independent of direc- 
tion, though it requires a higher field than reached in the 
current experiment to say for certainty that saturation 
has been achieved for B||[110]. We have no definitive 
explanation for the lack of anisotropy in the measured 
saturation moment, perhaps the influence of the dipolar 
interaction, which is not taken into account by Harris et 
al. ||, plays a role. We should also point out that the 
calculations of Harris et al. || are at T/J = 0.1 ~ 0.3if, 
which is much lower than our lowest temperature, but 
the qualitative agreement is good since both tempera- 
tures are below T* . This is the same conclusion of a 
different simulation that shows a single step like feature 
followed by saturation for B||[lll] and an identical sat- 
uration moment (though a much different M versus B 
curve) for a different field direction jl7]] . 

We have performed the first measurements of the mag- 
netization on Ho2Ti207 single crystals along various 
crystalline directions. The magnetization results can be 
explained by S = 1/2 spins with g = 20 which lock in 
along the [111] direction at T* » 1.97 K, and the results 
of applying the field along the different crystalline axes 
are explained qualitatively by a Monte Carlo calculation 
of the spin-ice model with nearest neighbor interactions 
only ||. In particular, the observation of a plateau in 
the magnetization as a function of applied magnetic field 
for 5 1 1 [111] at T < T* as predicted [§, has been ob- 
served for the first time. Contrary to predicted behavior 
11, it was found that the saturation magnetic moment 
of ~ 5.9 /ib/Ho is nearly independent of the orientation 
of the crystalline axes relative to the applied field. The 
magnetization and heat capacity measurements are both 
consistent with an interaction strength of J w 7.2 — 7.4 
K between the Ho magnetic moments. 

In conclusion, above T* the spins have a preferred 
orientation along the [111] direction due to single ion 
anisotropy revealed by the reduced saturated moment 
in the system. From magnetization measurements be- 
low 2 K the moments on a single tetrahedron have a 
specific two-spins-in, two-spins-out configuration and the 
moments cannot overcome an energy barrier to another 
configuration. This result is consistent with neutron scat- 
tering work which finds the Ho spins in a two-spins-in, 
two-spins-out configuration and recent specific heat re- 
sults g. 
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FIG. 2. The inverse static susceptibility in an applied 
field of 1 kOe as a function of temperature. The line is a 
Curie- Weiss fit as described in text. The inset shows the 
magnetic entropy AS m as a function of temperature T for 
Ho2Ti2 07 at T and 1 T. The solid line corresponds to the 
expected values of the nearest neighbor spin-ice model as dis- 
cussed in the text. The dashed lines represent the recovery of 
(l/2)i?ln(3/2) of the expected 7?ln2 entropy as discussed in 
text. 



FIG. 3. Magnetization as a function of applied magnetic 
field along the listed crystalline directions for Ho2Ti2 07 at 
various temperatures. A plateau is clearly observed in the 
data for B||[lll] as the temperature is lowered below T* ~ 2 
K. 



FIG. 1. AC susceptibility times temperature versus inverse 
temperature on a Ho2Ti2 07 single crystal for magnetic field 
applied along the [111] axis. The horizontal line is the ex- 
pected value for g — 20. The other line is a fit described in 
the text. 
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Cornelius and Gardner, Fig. 1 
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